We  introduce  a  version  of  the  asymptotic  expansions  for  Bessel  functions  Ju(z),  Yv(z) 
that  is  valid  whenever  \z\  >  is  (which  is  deep  in  the  Fresnel  regime),  as  opposed  to 
the  standard  expansions  that  are  applicable  only  in  the  Fraunhofer  regime  (i.e.  when 
1 2 1  >  is2).  As  expected,  in  the  Fraunhofer  regime  our  asymptotics  reduce  to  the  clas¬ 
sical  ones.  The  approach  is  based  on  the  observation  that  Bessel’s  equation  admits 
a  non-oscillatory  phase  function,  and  uses  classical  formulas  to  obtain  an  asymptotic 
expansion  for  this  function;  this  in  turn  leads  to  both  an  analytical  tool  and  a  nu¬ 
merical  scheme  for  the  efficient  evaluation  of  Ju(z),  Yv(z),  as  well  as  various  related 
quantities.  The  effectiveness  of  the  technique  is  demonstrated  via  several  numerical 
examples.  We  also  observe  that  the  procedure  admits  far-reaching  generalizations  to 
wide  classes  of  second  order  differential  equations,  to  be  reported  at  a  later  date. 
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On  the  asymptotics  of  Bessel  functions  in  the  Fresnel 

regime 


Given  a  differential  equation 


y"{z)  +  q(z)y{z)  =  0,  (1) 

a  sufficiently  smooth  a  :  C  — >  C  is  referred  to  as  a  phase  function  for  (1)  if  the  pair  of 
functions  u,  v  defined  by  the  formulae 


COs(q:(^)) 

(2) 

\a'(z )|1/2 

sin(a(,s)) 

(3) 

\a'(z)\1/2 

forms  a  basis  in  the  (two-dimensional)  space  of  solutions  of  (1). 

Phase  functions  arise  from  the  theory  of  global  transformations  of  ordinary  differential 
equations,  which  was  initiated  in  [5];  more  modern  discussions  can  be  found,  inter  alia,  in 
[2],  [6].  Despite  their  long  history,  phase  functions  possess  a  property  that  appears  to  have 
been  overlooked:  as  long  as  the  function  q  is  non-oscillatory,  the  equation  (1)  possesses  a 
non-oscillatory  phase  function.  This  observation  in  its  full  generality  is  somewhat  technical, 
and  will  be  reported  at  a  later  date;  in  this  short  note,  we  apply  it  to  the  case  of  Bessel’s 
equation 

Z2(p"(z)  +  zip'(z)  +  (z2  -  is2)<p(z )  =  0.  (4) 

To  the  authors’  knowledge,  the  results  in  this  paper  are  new  despite  the  long-standing 
interest  in  the  phase  functions  associated  with  Bessel’s  equation.  The  most  closely  related 
antecedents  of  this  work  appear  to  be  [3]  and  [9,  10,  7].  The  former  uses  Taylor  expansions 
of  a  non-oscillatory  phase  function  for  Bessel’s  equation  to  evaluate  Bessel  functions  of 
large  arguments;  and  the  latter  three  works  make  use  of  phase  functions  to  compute  zeros 
of  Bessel  functions  (as  well  as  certain  other  special  functions). 

1  Phase  functions  and  the  Kummer  Equation 

This  section  contains  a  summary  of  well-known  facts,  all  of  which  can  be  found  in  [2,  6]. 
We  say  that  a  C3  function  a(z)  is  a  phase  function  for  the  differential  equation  (1)  if  the 
functions  u,  v  defined  by  the  formulae 

u(z)  =  \J u2(z)  +  v2{z)  cos(a(z)),  (5) 

v(z)  =  \Ju2(z)  +  v2(z)  sin(a(^))  (6) 
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are  solutions  of  (1).  All  phase  functions  for  (1)  are  solutions  of  the  third  order  nonlinear 
ordinary  differential  equation 


=  q(t) 


1  a'"(t)  3  (a"(t)\2 

2  a'(t)  4  \  a'(t)  ) 


(7) 


Likewise,  any  solution  of  (7)  yields  a  phase  function  for  the  differential  equation  (1);  we 
will  refer  to  (7)  as  Kummer’s  equation,  after  E.  E.  Kunimer  who  studied  it  in  [5].  Another 
simple  connection  between  the  solutions  u,v  of  (1)  and  the  corresponding  phase  function 
is  given  by  the  formula 

,,  v  w 

u2(z)  +  v2(zy 

with  W  the  (constant)  Wronskian  of  the  basis  {u,v}. 

Our  principal  interest  is  in  the  highly  oscillatory  case,  where  q(z)  is  of  the  form  7 2q(z) 
with  7  a  large  real-valued  constant  and  q  a  complex-valued  function,  so  that  q(z)  is  asymp¬ 
totically  of  the  order  y2.  As  a  consequence  of  the  Sturm  comparison  theorem  (see,  for 
example,  [4])  solutions  of  (1)  are  necessarily  highly  oscillatory  when  7  is  large.  But,  in- 
triguingly,  (5)  and  (6)  leave  open  the  possibility  that  a  non-oscillatory  phase  function  a(z) 
can  be  used  to  represent  highly  oscillatory  functions  u  and  v:  in  fact,  such  a  are  available 
whenever  q  is  non-oscillatory  (in  an  appropriate  sense).  In  this  note,  we  construct  such  a 
in  the  case  of  Bessel’s  equation. 


2  A  Phase  function  for  Bessel’s  equation 


The  Bessel  function  of  the  first  kind  of  order  v 


J,Az )  =  (  - 


£ 


(-P2) 


^  n\  r  (u  +  n  +  1) 

n= 0 


and  the  Bessel  function  of  the  second  kind  of  order  v 

Ju{z)  cos(z/7 t)  —  J-U(z) 


YJz  = 


Sin  U/7T 


form  a  basis  in  the  space  of  solutions  of  Bessel’s  equation  (4). 

Bessel  functions  are  among  the  most  studied  and  well-understood  of  special  functions. 
Even  the  standard  reference  books  such  as  [1],  [4]  contain  a  wealth  of  information.  Here, 
we  only  list  facts  necessary  for  the  specific  purposes  of  the  paper;  the  reader  is  referred 
to  [11],  [1],  [4]  (and  many  other  excellent  sources)  for  more  detailed  information. 

For  arbitrary  positive  real  v  and  large  values  of  |z|,  the  functions  ./„,  Yv  possess  asymp¬ 
totic  expansions  of  the  form 


Ju{z) 


7T 


4 


(9) 


3 


Y„(z) 


Z/7T 


sin  z 


7 r 


4 


k= 1  ~  / 


(10) 


where  for  each  k  =  1,  2,  •  •  • ,  P2fc,  are  polynomials  of  order  2k  (see,  for  example,  [1]  for 
exact  definitions  of  {P2k},  {Q2k})- 

Remark  2.1.  As  one  would  expect  from  the  form  of  expansions  (9),  (10),  they  provide  no 
useful  approximations  to  Jv(z),  Yv(z)  until  \z\  is  greater  than  roughly  v2 .  The  regime  where 
v  <  \z\  <  v2  is  known  as  Fresnel  regime,  while  the  regime  \z\>  v2  is  referred  to  as  Fraun¬ 
hofer  regime.  In  other  words,  in  the  Fraunhofer  regime  the  classical  asymptotics  (9),  (10) 
become  useful. 

The  transformation  i/)(z)  =  z1^2(p(z)  brings  (4)  into  the  standard  form 


r(z)  + 


if(z)  =  0. 


(11) 


Two  standard  solutions  of  (11)  are  {y/zJu(z),  y/zYv(z)};  in  this  case,  (8)  becomes 

2  1 


a„(z  = 


KZ  J2{z)  +Y?(z) 

Clearly,  (12)  determines  the  phase  function  a  up  to  a  constant;  choosing 


7T 


atv(z)  =  --  +  /  a'u(u)du, 


Jo 


we  obtain  expressions 


Ju(z)  =  Mv(z)  cos(av(2)) 


(12) 


(13) 

(14) 


Yv(z)  =  Mv(z)sm(av(z)),  (15) 

with  Mv(z)  defined  by  the  formula 

Mv(z)  =  y/j*{z)  +  Y*(z).  (16) 


In  a  remarkable  coincidence,  there  exists  a  simple  integral  expression  for  Mv,  valid  for 
all  z  such  that  arg(z)  <  n.  Specifically, 


o  roc 

(Mu(z))2  =  J2(z)  +Y2(z)  =  —  /  K0(2zsmh.(t))  cosh(2i/f)  dt  (IT) 

71  Jo 

(see,  for  example,  [4],  Section  6.664);  even  a  cursory  examination  of  (17)  shows  that  for  z 
on  the  real  axis,  Mu  is  a  non-oscillatory  function  of  z. 

The  approximation 


Ml(z) 


2  ^  T  (n  +  1)  T  (y  +  1  +  n)  1 
nz  n\Jn  T  (u  +  \  —  n)  z2n  ’ 

n= 0  v  \  z  / 


(18) 
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can  be  found  (for  example)  in  Section  13.75  of  [11];  it  is  asymptotic  in  and  its  first 
several  terms  are 


M2v{z)~—  1  +  - 


7TZ 


l/i-l  1  3(/i-l)(/i-9)  1  3  5(/i-  l)(/i  —  9)(/x  —  25) 

V  +  2  '  4  (2z)4  +  2  '  4  '  6  pi)5  + 

(19) 

with  //  =  4z/2.  When  the  expansion  (18)  is  truncated  after  k  terms,  the  error  of  the  resulting 
approximation  is  bounded  by  the  absolute  value  of  the  next  term,  as  long  as  z  is  real  and 
k  >  is  (see  [11]).  Computationally,  it  is  often  convenient  to  rewrite  (18)  in  the  form 


M2M 


2 

7 TZ 


(20) 


with  to  —  1,  and 


f n~  (2n-  1)2\  2n  —  1 
V  4  )  2  n 


(21) 


Observation  2.1.  While  the  expansion  (19)  has  been  known  for  almost  a  century,  one 
of  its  implications  does  not  appear  to  be  widely  understood.  Specifically,  (19)  provides  an 
effective  approximation  to  Mv  whenever  \z\  >  is;  combined  with  (If),  (15),  it  results  in 
asymptotic  approximation  to  JU,YU  in  the  Fresnel  regime,  i.e.  when  \z\  >  is  -  as  opposed  to 
the  expansions  (9),  (10),  that  are  only  valid  in  the  Fraunhofer  regime  -  i.e.  when  \z\  >  is2 . 
This  (rather  elementary)  observation  has  obvious  implications,  inter  alia,  for  the  numerical 
evaluation  of  Bessel  functions  of  high  order. 

Applying  Lemma  5.1  to  the  expansion  (20)  and  using  (12),  (16),  we  obtain 


a„(z)  = 


2  1 

7 TZ  Mf{z) 


y2n 


(22) 


71=1 


with  so  =  1,  and  si,  S2,  S3, . . .  defined  by  the  formula 


71—1 

Sn  =  —  I  Sn  +  ^  ^  t 
3= 1 


JSn-J  I  1 


(23) 


the  first  several  terms  in  (23)  are 

.  .  .  ,  fi-  1  p2  —  26/i  +  25  /a3-  115/i2  +  1187/r  -  1073 

d  ( z )  X  — - — - —  - _ I- 

uK  ’  8^2  128z4  1024^6 

Obviously,  the  indefinite  integral  of  (22)  is 


(24) 


a.v(z)  ~  c + z  - 


— [  (2 n  —  1)  z2n  1  ’ 


(25) 
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with  C  to  be  determined;  the  first  several  terms  in  (25)  are 

,  N  ^  u  —  1  u2  —  26n  +  25  u3  -  115/r2  +  1187/i  -  1073 
a^)~C  +  Z+V+  384,3  +  ^ - +  "'•  <26> 

In  order  to  find  the  value  of  C  in  (25),  (26),  we  observe  that  for  sufficiently  large  \z\,  the 
expansion  (26)  becomes 

av(z)  ~C  +  z.  (27) 

Substituting  (27)  into  (14),  (15),  we  have 

Jv{z)  ~  Mv{z)  cos (C  +  z)f  (28) 

Yv(z)  ~  M„(z)  sin (C  +  z).  (29) 

Clearly,  the  approximations  (28),  (29)  must  be  compatible  with  (9),  (10),  which  yields 

(30) 

Finally,  substituting  (30)  into  (26),  we  end  up  with 


_  Z/7T  7T 

(  “  2  4 


VTT  TT 

0‘u(z)^-—--+Z~^ 


— '  (2n  —  1)  z2n_1  ’ 


(31) 


or 


ay(^) 


7T  /i  -  1  yU2  —  26/r  +  25  yU3  -  115/i2  +  1187/r  -  1073 
4  +  ~ +  8^  +  384^3  +  5120^5  + 


(32) 


with  fi  =  4v2.  We  note  that  the  asymptotic  expansion  (31)  is  only  useful  because  the 
particular  phase  function  av(z)  it  represents  is  non-oscillatory:  Figure  3  shows  plots  of  the 
derivative  of  the  phase  function  associated  with  {^Ju(z),  \fzYv(z)}  and  also  the  derivative 
of  the  phase  function  associated  with  another  choice  of  basis,  {2 yfzJv{z),  ^Yu(z)}. 


3  Numerical  Experiments 


In  this  section  we  present  the  result  of  several  numerical  experiments  conducted  to  verify 
the  scheme  of  this  paper.  The  code  for  these  experiments  was  written  in  Fortran  77  and 
compiled  using  the  Intel  Fortran  compiler  version  12.0.  Experiments  were  conducted  on  a 
laptop  equipped  with  an  Intel  Core  i7-2620M  processor  running  at  2.70  GHz  and  8  GB  of 
RAM.  Machine  zero  was  £$  =  2.22044604925031  x  10“16. 

We  constructed  asymptotic  approximations  to  the  functions  ./„,  Yv  via  the  formulae  (21),  (23),  (31). 
In  order  to  avoid  exceeding  the  machine  exponent,  we  altered  the  procedure  slightly,  so 
that  the  coefficients  tk,  Sk  are  never  computed  by  themselves:  only  the  ratios 


tk  Sk 
~2 k  ’  ~2 k 


(33) 
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are  calculated  via  obvious  modifications  of  (21),  (23). 

As  with  all  procedures  relying  on  asymptotic  expansions,  it  is  not  always  possible  to 
achieve  a  desired  accuracy.  Indeed,  the  magnitudes  of  the  terms  of  the  expansion  reach  a 
certain  minimum  and  then  proceed  to  increase.  And,  of  course,  truncating  the  expansions 
when  the  terms  become  small  does  not  necessarily  ensure  the  accuracy  of  the  approximation 
(see  [8]  for  numerous  examples  of  possible  pathologies).  It  is  shown  in  Section  13.75  of  [11] 
that  if  v  is  real,  z  is  positive  and  n>  v  —  1/2,  then  the  remainder  resulting  from  the  first 
n  terms  of  expansion  (19)  is  smaller  in  magnitude  than  the  (n  +  l)st  term.  However,  the 
authors  are  not  aware  of  any  error  bounds  for  (31)  (or  for  (19)  )  in  the  general  case. 

Moreover,  the  value  of  the  phase  function  av(z)  is  proportional  to  the  argument  z  and 
values  of  Ju(z)  and  Yu(z)  are  obtained  in  part  by  evaluating  the  sine  and  cosine  of  av(z). 
This  imposes  limitations  on  the  accuracy  of  the  obtained  approximations  when  z  is  large 
due  to  the  well-known  difficulties  in  evaluating  periodic  functions  of  large  arguments. 

3.1  Comparison  with  Mathematica 

In  these  first  experiments,  we  applied  the  procedure  of  this  paper  to  the  evaluation  of  the 
Bessel  functions  Jv{z)  and  Yv(z)  at  various  orders  and  arguments.  The  resulting  values 
were  compared  with  those  produced  by  version  9.0.0  of  Wolfram’s  Mathematica  package; 
30  digit  precision  was  requested  from  Mathematica.  Table  2  reports  the  results.  There, 
the  number  of  terms  used  in  the  expansions  of  the  modulus  and  phase  functions  and  the 
relative  errors  in  the  obtained  values  of  Jv(z)  and  Yv(z)  are  reported. 

3.2  Bessel  functions  of  large  order. 

In  this  experiment,  we  approximate  the  values  of  Bessel  functions  of  very  large  orders  and 
arguments.  Comparison  with  other  approaches  is  difficult  for  such  large  orders;  for  instance, 
Mathematical  Bessel  function  routines  are  prohibitively  slow  in  this  regime.  We  settled 
for  running  our  procedure  twice,  once  using  double  precision  arithmetic  and  once  using 
extended  precision  (Fortran  REAL*16)  arithmetic  in  order  to  produce  reference  values  for 
comparison. 

The  first  row  of  each  entry  in  Table  3  reports  the  relative  error  in  the  approximations 
of  Ju(z)  and  the  second  row  gives  the  relative  error  in  the  approximation  of  Yl/(z).  Table  1 
gives  the  number  of  terms  in  the  expansion  of  the  modulus  function  and  the  number  of 
terms  in  the  expansions  of  the  modulus  and  phase  function  used  to  evaluate  Jv{z)  and 
Yv{z).  These  values  depended  only  on  the  ratio  of  z  to  v  and  not  on  the  value  of  v. 

3.3  Failure  for  small  orders. 

In  this  experiment,  we  considered  the  performance  of  the  procedure  of  this  paper  for  rela¬ 
tively  small  values  of  \v\.  We  evaluated  J„(  10v)  at  a  series  of  values  of  u  between  0  and  5. 
A  plot  of  the  base-10  logarithm  of  the  relative  error  in  J„(10//)  is  shown  in  Figure  2.  Errors 
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were  estimated  via  comparison  with  Wolfram’s  Mathematica  package;  30  digit  precision 
was  requested  from  Mathematica. 

3.4  Failure  as  arg(z)  approaches  7 r. 

In  this  experiment,  we  computed  the  values  of  Yio(exp(id))  as  6  approaches  7 r.  The  ob¬ 
tained  values  were  compared  to  those  reported  by  Wolfram’s  Mathematica  package;  30  digit 
precision  was  once  again  requested  from  Mathematica.  A  plot  of  the  base-10  logarithm  of 
the  relative  error  in  Wo  (exp  (id))  as  a  function  of  6  is  shown  in  Figure  1. 

4  Conclusions 

We  have  shown  that  the  Bessel  functions  Ju{z)  and  Yu(z)  can  be  efficiently  evaluated  when 
v  is  large  and  \z\  >  \v\.  This  was  achieved  by  representing  Bessel  functions  in  terms  of 
a  non-oscillatory  phase  function  for  which  (conveniently  enough)  a  well-known  asymptotic 
expression  is  available. 

The  observation  underlying  the  scheme  of  this  paper  —  namely,  the  existence  of  a 
non-oscillatory  phase  function  —  is  not  a  peculiarity  of  Bessel’s  equation.  The  solutions 
of  a  large  class  of  second  order  linear  differential  equations  can  be  approximated  to  high 
accuracy  via  non-oscillatory  phase  functions,  a  development  the  authors  will  report  at  a 
later  date. 


5  Appendix 


Here,  we  formulate  a  lemma  used  in  Section  2;  its  proof  is  an  exercise  in  elementary  calculus, 
and  can  be  found,  for  example,  in  [8]. 

LEMMA  5.1.  Suppose  that 


f(z)  ~ 


cq  a-2  03 

lH - 1 - 1 — -  4- 

z2  z4  z6 


(34) 


is  an  asymptotic  expansion  for  f  :  C  — >•  C ,  with  cq,  cq,  <23, . . .  a  sequence  of  complex  num¬ 
bers.  Then  the  asymptotic  expansion  of  1  / f  is 


1 


-  rs_/ 

m 


i+V 

z 2 


bo  b:> 

4  +  4  +  --- 

z4  z6 


(35) 


where  b\  =  —  eq,  and  the  rest  of  the  coefficients  tq,  63,  64, . . .  are  given  by  the  formula 


n—  1 

—  ajbn_j. 
l=i 

8 


bn  & n 


(36) 
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zjv  Modulus  terms  Phase  terms 


1.1 

162 

133 

2 

024 

022 

10 

008 

008 

107T 

006 

006 

Table  1:  Bessel  functions  of  large  order.  The  number  of  terms  in  the  expansions  of  the 
modulus  and  phase  functions  used  in  the  computation  of  Bessel  functions  of  large  orders. 


Figure  1:  Failure  as  arg(z)  approaches  i r.  The  base-10  logarithm  of  the  relative  error 
in  the  approximation  of  y10(100exp(f61))  as  6  approaches  tt. 
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Order  v 

Figure  2:  Failure  for  small  orders.  The  base-10  logarithm  of  the  relative  error  in  the 
approximation  of  Ju(  10u)  plotted  as  a  function  of  v. 


Figure  3:  On  the  left,  a  plot  of  the  derivative  of  a  phase  function  associated  with  the 
basis  of  solutions  {\/zJ2o(z),  y/zY2o(z)}  for  (11);  on  the  right,  a  plot  of  the  derivative  of  a 
phase  function  associated  with  the  basis  {2^/zJ2q{z)1  ^Y20(z)}.  Both  are  over  the  interval 
[0,200], 
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Modulus  Phase  Relative  error  Relative  error 


V 

2 

terms 

terms 

in 

Jv{z) 

in 

Y„(z) 

50 

1.1 

v\ 

049 

091 

1.45 

X 

10- 

-14 

3.25 

X 

10- 

-14 

2 

v\ 

023 

022 

3.73 

X 

io- 

-14 

9.80 

X 

10" 

-15 

10 

v\ 

008 

008 

1.77 

X 

io- 

-14 

9.38 

X 

10- 

-15 

100| 

v\ 

005 

005 

4.53 

X 

io- 

-13 

7.23 

X 

io- 

-14 

50  -  lOz 

l.ll 

v\ 

075 

098 

4.85 

X 

io- 

-15 

4.81 

X 

10" 

-15 

2 

v\ 

023 

022 

6.04 

X 

io- 

-15 

6.12 

X 

10“ 

-15 

10 

v\ 

008 

008 

9.79 

X 

io- 

-14 

9.77 

X 

10" 

-14 

100| 

v\ 

005 

005 

6.31 

X 

io- 

-13 

6.31 

X 

10“ 

-13 

100  +  20z 

IT 

v\ 

080 

117 

1.79 

X 

io- 

-14 

1.79 

X 

10“ 

-14 

2 

v\ 

024 

021 

3.19 

X 

io- 

-14 

3.19 

X 

10" 

-14 

10 

v\ 

008 

008 

1.46 

X 

io- 

-13 

1.46 

X 

10“ 

-13 

100| 

v\ 

005 

005 

5.43 

X 

io- 

-13 

5.43 

X 

10" 

-13 

10000 

IT 

v\ 

162 

133 

1.72 

X 

io- 

-11 

4.63 

X 

10“ 

-12 

2 

v\ 

024 

022 

1.18 

X 

io- 

-12 

6.69 

X 

10“ 

-13 

10 

v\ 

008 

008 

8.96 

X 

io- 

-13 

1.28 

X 

10" 

-10 

100| 

v\ 

005 

005 

1.10 

X 

io- 

-10 

2.96 

X 

10“ 

-12 

100000 

IT 

v\ 

162 

133 

5.18 

X 

io- 

-13 

7.16 

X 

10" 

-13 

2 

v\ 

024 

022 

4.82 

X 

io- 

-11 

2.34 

X 

10“ 

-11 

10 

v\ 

008 

008 

1.57 

X 

io- 

-09 

3.93 

X 

10“ 

-12 

100| 

v\ 

005 

005 

1.92 

X 

10- 

-10 

2.41 

X 

10" 

-10 

Table  2:  Comparison  with  Mathematica.  The  results  of  the  experiments  of  Section  3.1. 
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v  =  106  V  =  109  v  =  1012  1/  =  1015  1/  =  1018 


1  It. 

1.96 

X 

10- 

-li 

4.48 

X 

10- 

-07 

3.74 

X 

10- 

-07 

3.26 

X 

10- 

-08 

5.78 

X 

10- 

-06 

1.1  V 

1.18 

X 

io- 

-10 

2.81 

X 

io- 

-07 

1.70 

X 

io- 

-08 

1.57 

X 

io- 

-06 

5.87 

X 

10" 

-08 

1.04 

X 

10" 

-10 

2.15 

X 

io- 

-07 

1.01 

X 

io- 

-07 

1.25 

X 

io- 

-07 

4.08 

X 

10" 

-07 

4.66 

X 

10" 

-11 

6.94 

X 

io- 

-07 

1.37 

X 

io- 

-07 

2.96 

X 

io- 

-07 

7.73 

X 

10" 

-07 

5.25 

X 

10" 

-09 

3.80 

X 

io- 

-07 

3.20 

X 

io- 

-05 

1.09 

X 

io- 

-05 

1.04 

X 

10" 

-06 

1UZ/ 

4.53 

X 

10" 

-10 

5.30 

X 

io- 

-07 

4.33 

X 

io- 

-09 

6.83 

X 

io- 

-08 

4.00 

X 

10" 

-06 

3.23 

X 

10" 

-10 

4.58 

X 

io- 

-06 

7.33 

X 

io- 

-06 

3.39 

X 

io- 

-06 

4.12 

X 

10" 

-06 

1U7 TV 

3.80 

X 

10- 

-10 

1.01 

X 

10" 

-05 

2.10 

X 

10“ 

-07 

3.68 

X 

10" 

-06 

3.05 

X 

10- 

-06 

Table  3:  Bessel  functions  of  large  order.  The  relative  errors  in  the  obtained  approxi¬ 
mations  of  Bessel  functions  of  large  orders:  the  round-off  due  to  large  values  of  \z\  sharply 
limits  the  accuracy  of  results  obtainable  in  double  precision 
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